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' Abstract. We estimated the magnetic field strength of compact stars in X-ray binaries 
• • Sco X-1 and Cyg X-1, via various methods of determination of magnetic fields. For Sco 
■"^ • X-1 we used three independent methods. One of them is based on the correct account of 
the Faraday rotation of polarization plane in the process of electron scattering of X-rays 



from accreting neutron stars. Numerical calculations are made with use of first X-rays po- 
larimetric data presented by Long et al. (1979). Other original method of determing the 
magnetic field developed by Titarchuk at al. (2001) is based on observed quasi-periodic 
oscillations (QPO) frequencies in X-ray binaries that can be considered as magnetoa- 
coustic oscillations of boundary layer near a neutron star. The optical polarimetric data 
obtained in 70-th have been also used for estimation of magnetic field of the neutron star 
in Sco X-1 and of nearest environment around the black hole in Cyg X-1. 
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1. Introduction 

For X-ray sources of type Sco X-1 the Compton scattering by electrons is very important 
for generation of the radiation spectrum. The Compton scattering produces the Unear 
polarization of radiation from such type objects. The polarization of X-ray sources have 
been first considered by Rees (1975) and independently by Lightman & Shapiro (1975). 
They calculated the spectrum of polarized radiation from an accretion disk. 

But the spectrum of polarized radiation is drastically changed due to the Faraday 
rotation of the polarization plane in a magnetic field. Gnedin & Silant'ev (1980, 1997) 
(see also the book by Dolginov et al. 1995) have calculated the spectrum of polarized 
radiation for the various astrophysical objects with magnetic fields taking into account 
the Faraday rotation effect and suggested the new method of determination of stellar 
magnetic field with the spectrum of broad band linearly polarized radiation (see, for 
example, Gnedin & Silant'ev 1984; Gnedin et al. 1988). The application of their method 
to accretion disks around supermassive black holes (AGNs, QSOs) has been recently 
developed by Gnedin & Silant'ev (2002). 

The first rough estimation of a magnetic field of Sco X-1, based on the first X-ray 
polarimetric observation made by Long et al. (1979), was done by Gnedin & Silant'ev 
(1980). 

Recently Titarehuk et al. (2001) presented the detail model of quasi-periodic oscilla- 
tions (QPO) of Sco X-1 using a large set of Rossi X-ray Timing Explorer (RXTE) data 
for this object. The QPOs have been observed not only into kHz range, but also at 6 
Hz. There is an almost linear correlation between Sco X-l's kHz and 6 Hz QPOs. They 
identify the 6 Hz frequencies as acoustic oscillations of a Keplerian disk around the 
neutron star (NS) that is formed after radiation pressure near the Eddington accretion 
rate destroys the disk. They showed that the disk is placed near one NS radius from the 
siirfaee of NS. It is the remarkable fact that they could estimate the Sco X-1 magnetic 
field. It appears to be 0.7 x lO^G at about one NS radius above the NS surface. 

Let us remind also that still Gnedin & Shulov (1971) and Nikulin et al. (1971) claimed 
some evidence of circularly polarized light of Sco X-1 in optical spectral range. Though 
the followed measurement (see Illing & Martin 1972) did not confirm the existence of 
circular polarization for Sco X-1, Kemp et al. (1972) claimed the possible existence of 
the strong nonstationary circular polarization state. They called this state as a "flaring 
polarization" . Nevertheless the situation with determination of magnetic field of Sco X-1 
was unclear till now. The last results by Titarchuk et al. (2001), allow to reconsider this 
situation and to compare the results of various methods of magnetic field estimations. 
We start with the detail description of presented methods of magnetic field estimations 
of compact stars emitting X-rays. 
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2. Linear Polarization Spectrum of X-ray Radiation from Magnetized Disks 
and Envelopes Around a Compact Star 

The X-ray radiation of a compact star (neutron star, black hole), scattered on electrons 
of the envelope of an accreting matter, acquires linear polarization. It is well known, 
however, that integral polarization from a spherical envelope is equal to zero due to 
total compensation of directions of electric vectors from various areas of this envelope. If 
magnetic field is absent, the nonzero net polarization requires nonspherical shape of the 
scattering electron envelope. If there is the presence of a magnetic field in the circumstellar 
electron envelope, the situation is radically changed. The scattered radiation undergoes 
the Faraday rotation of its polarization plane due to the propagation in a magnetized 
plasma. 



The angle of the Faraday rotation x is determined by the expression (see Gnedin & 
Silant'ev 1980, 1997) 

X = 1;6tt cose, 



where tt is the Thomson optical thickness of the envelope, 9 is the angle between the 
directions of a magnetic field B and radiation propagation n, = e^/mgC^ is classic 
electron radius, lub = eB/meC is the cyclotron frequency. 

If the magnetic field increases, the angle x increases too and partly polarized scattered 
radiation begins to undergo the Faraday rotation. The rotation angles x ^re different for 
waves scattered in different volumes along the line of sight n, and, as a result, the total 
radiation from all volume elements will be depolarized (see Fig. 10 from the review by 
Gnedin & Silant'ev 1997). 

The spectrum of polarized radiation from a spherically symmetric envelope has been 
calculated by Gnedin & Silant'ev (1980). For a spherical envelope with the dipole mag- 
netic field only the magnetic equator volumes don't produce the Faraday rotation of 
polarization plane because of for these volumes cos9 w and x ~ 0. As a result, only the 
radiation scattered in equatorial volumes acquires net polarization in the plane (nM), 
where n and M are the directions of line of sight and dipole magnetic field, correspond- 
ingly. Namely this radiation alone gives rise to non-zero integral polarization from the 
magnetized envelope. The shape of the spectrum of the integral linear polarization has 
a bell-shaped form. Indeed, for very short wavelengths the Faraday rotation angles tend 
to zero and the integral polarization disappears due to spherical symmetry of the enve- 
lope. For very long wavelengths the angles of Faraday's rotation increase and the near 
equatorial region with x ~ decreases what diminishes the integral polarization. So, the 
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maximum integral polarization corresponds to such magnetic field strengths and wave- 
lengths when the mean rotation angle x ~ 1- The process of strong depolarization begins 
with the the values of magnetic field and wavelengths when 6tt > 1 (see Eq. (1)). 

For disk-like electron envelope without the true absorption with tt ^ 1 we have for 
outgoing radiation (see Silant'ev 2002): 

Pi(n, B) ^ 0.914 ^ (2) 

where = cosi?, and is the angle between the line of sight n and the normal to the 
disk surface. The terms without the square root describe roughly the polarization degree 
P;(n) of radiation in the absence a magnetic field. So, the formula can be written as 

P,(n,B)^ /iTa^"^ 2, - (3) 
Vl + cos^ f) 

Later we shall use this formula for the estimation of magnetic field in the environment 
of the source Sco X-1. 

For the spherically symmetric optically thin envelope around the neutron star with 
the dipole magnetic field the degree of linear polarization of X-rays scattered into this 
envelope can be estimated by the expression 

Pi(n, B) = TTf{Ro/Rs; SsTo; 9) sin^ 0. (4) 

Here Ro/Rs is the ratio of the envelope radius to the neutron star radius, Sg is the 
depolarization parameter (1) that corresponds to the stellar equatorial magnetic field 
strength. The function / is tabulated by Gnedin & Silant'ev (1980, 1984). 

In the region of the depolarization when 6sTt > 1 the asymptotic expression of the 
function / takes a place: 

/ « {SsTrr^ ~ {u/ojo)^^. (5) 

Here n is the exponent of radial dependence of the electron density in the envelope: 
A^e ~ {n 7^ 0). The value too = (3a;Bc/2re)^/^. If, for example, n = 2, the function 

/ drops with frequency decrease as / w^/^. For 5sTt <C 1 (large frequencies) the 
spectrum of polarization degree has universal dependence: Pi ~ W*. 

3. Magnetic Field Estimation from Magnetoacoustic Oscillations in Neutron 
Star Binaries 

The original method to determine the magnetic field around neutron stars based on 
observed kHz and viscous quasi-periodic oscillations (QPO) frequencies has recently been 
proposed by Titarchuk et al. (2001). They have analyzed magnetoacoustic wave formation 
on the layer between a neutron star surface and the inner edge of a Keplerian disk and 



Yu. Gnedin, N. Silant'ev & L. Titarchuk: Determination of magnetic field 



derived formulas for the magnetoacoustic wave frequencies for different regimes of radial 
transition layer oscillations. As a result they demonstrated that one can use the QPO 
as a new kind of probe to determine the magnetic field strengths of neutron stars in the 
binaries. 

For the two extreme (acoustic and magnetic) cases Titarchuk et al. (2001) have been 
obtained an approximate formula for the magnetoacoustic frequency coma' 

. 

(^)V| (^)^(^)Vi(w)" 



+ 



(6) 



Here Vs is the sound velocity, Va is the Alfven velocity, rout and ri„ are the outer and 
the inner radii of a transit layer. The index a is related to the multipole magnetic field 
through the expression of V| r~", so that a = 6, 8, 10 are for the dipole, quadrupole, 
and octupole, respectively. 

The parameters /3s and Pm are determined by the transcendental equation (see 
Titarchuk et al. 2001): 
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where /3 = Zout - Zin and r] = Zout/zin = (?"o«t/ri„)("+^)/^, Zout and Zi„ are z-axis for the 
transit layer. 

Titarchuk et al. (2001) estimated approximately the values of (is and /3m: 



TT 2 

2+n 



/3m 

for a = 6 and 



/3s « {1.5/ [l + 1.57y/(ry-l)2]} 



(8) 



(9) 



for a = (the pure acoustic case). 

Substituting the sound velocity Vs = {kT/mpY/^ and the Alfven velocity Vm = 
-B (rout )/(47rp)^/^ one can estimate the magnetic field strength at the outer boundary of 
the transition layer. Then using the multipole law for the magnetic field one can also 
estimate the magnetic field strength at the neutron star surface. 



4. Magnetic Field Estimation from Circular Polarization of Optical Light of 
a Binary System 

The broad-band polarimctry of the optical continuum radiation is the effective direct 
method of the magnetic field estimation. The electromagnetic wave which is incident 
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onto plasma with a magnetic field produces oscillations of electron velocity. As a result 
an additional Lorentz force appears: 

F = f(^,Exb), (10) 

which depends on the ratio of the cyclotron frequency to radiation one and on the angle 
between the directions of wave's electric vector E and the magnetic field B. As a result of 
the action of this Lorentz force magnetized plasma acquires dichroism and birefringence 
properties, becoming similar to any anisotropic medium. In a magnetized plasma two 
types of electromagnetic waves (normal modes or normal waves) with different types of 
elliptical polarization should be propagated. These normal waves are usually called or- 
dinary (O.W.) and extraordinary (E.O.W.) with their intrinsic refraction indices, phase 
velocities and polarization. O.W. behaves as an usual electromagnetic wave in a plasma 
without a magnetic field. For E.O.W. transport coefficients for various emission (mag- 
netobremsstrahlung) and scattering (electron scattering) processes have resonance at 
cyclotron frequency wb (in detail see Dolginov et al. 1995). 

In the case when the radiation frequency w is much larger than the cyclotron frequency 
ujb, the radiation is predominantly circularly polarized: 

Pv = 2(wb/w) cosO; Pi ~ P^ ~ (wb/w)^. (11) 

One of the most important case is the cyclotron resonance: ui ojb (remember that 
ljb/oj 0.93- 10~^A(/xm)-B(G) ). For the optical range this case corresponds to magnetic 
field strengths B ~ 10^-j-lO^G. At cyclotron resonance the plasma radiation is completely 
polarized: 

„ sin^ 9^2 cos ^ „9 „9 , ^ 

1 + cos^ 6 1 + cos^ 6 

In our case of the close binary system with a magnetized neutron star the optical 
radiation of the accreting disk or plasma envelope around a neutron star can be circularly 
polarized. As a result the estimation of the neutron star magnetic field strength can be 
obtained from the measured optical circular polarization. 



5. Sco X-1: Estimation of Magnetic Field Strength 

The X-ray polarization of Sco X-1 was searched by Long et al. (1979). Multiple observa- 
tions of Sco X-1 were carried out with the Bragg crystal polarimcter aboard 0S08. An 
unambiguous detection of the time- averaged polarization was obtained for Sco X-1 which 
has a polarization Pi = (0.39 ± 0.20)% at 2.6 keV and Pi = (1.31 ± 0.40)% at 5.2 keV. 
These data give a signature of the polarization fall at energy values below 5.2 keV. If one 
suggests that this fall is due to the Faraday depolarization effect one can estimate the 
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magnetic field strength at an environment of the neutron star in Sco X-1 binary system. 
The ratio Pi{2.6keV)/Pii5.2keV) w 0.3. 

Using Eq. (3), we obtain for the 5.2 keV case ^cos6» ~ y(Pi(n)/1.31)2 - 1. The 
choice of Pi depends on of the angle of inclination of a disk. It is clear that we are to 
take Pi > 1.31%. As a result we have the estimation: 

B > 2.3 • 10^\/(Pi/1.31)2 - IG. (13) 

This estimation takes place for the model of optically thick inclined disk. If we choose 
Pi ~ 2, we obtain B > 2.6 ■ 10^ G. Analogous estimation for the 2.6 keV case gives 
B > 2.8 • 10*" G, i.e. in range of observational errors both estimations coincide one with 
another. It seems the polarized radiation forms at the inner part of accretion disk where 
the temperature is higher. 

Another estimation follows from the model of spherical electron envelope in dipole 
magnetic field. Considering that the polarization 1.31% corresponds to the position at 
the polarization spectrum curve which is near maximum, i.e. SgT/rf ~ 1, we obtain 
„3 

Bs > 2.3 • 10^ — G. (14) 

T 

Here, rj — Rq/Rs is the ratio of the envelope radius to the neutron star radius and r 
is the Thomson thickness of the envelope. Taking ~ 2 and r ~ 0.25, we obtain from 
(14) the estimatiom Eg ~ 7 ■ 10^ G, or the value B ~ 10® G in the region of an envelope. 
The estimations (14) considerably greater than those from (13). In reality this difference 
is more considerable because the condition Pi{5.2keV)/ Pi{2.6keV) ~ 3.3 can be carried 
out only if Sg > rj^ /t. 

Another way of estimation of Sco X-1 magnetic field is to use observed kHz and 
viscous QPO frequencies. This method was developed by Titarchuk et al. (2001). They 
used the best-fit parameters of Sco X-1 transition layer and determined the magnetic 
field strength Btl = (1-0 ± 0.05) • 10*^ G and Rtl = 2.12 • 10*^ cm in the transition layer. 
An extrapolation of the magnetic field from Rtl towards the neutron star radius gives 
us = 8 • 10^ G and = 3.3 x 10^ G for the dipole and octupole fields, respectively. 
These values appear close to the X-ray polarimetric estimation of neutron star magnetic 
field given by (13). 

Third way of estimation of Sco X-1 magnetic field is connected with observations of 
circular polarization of optical light in Sco X-1 binary system. 

Gnedin and Shulov (1971) and Nikulin et al. (1971) claimed occasional appearance of 
appreciable and variable circular polarization of optical light of Sco X-1 (see, also, Severny 
& Kuvshinov, 1975). Though lUing & Martin (1972) didn't confirm these results, Kemp 
et al. (1972) claimed also evidence of the occasional appearance of circular polarization 
of Sco X-1 that was named by them as "fiaring polarization". 
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The net circular polarization of Sco X-1 was estimated with Eq. (11). Gnedin & Shulov 
(1971) made the observations of circular polarization in the yellow-red light that corre- 
sponds the effective wavelength Ag// = 0.64/im. For their case estimation (11) give the 
following value of the magnetic field strength: 



This value corresponds quite well to the magnetic field magnitude derived by Titarehuk 
et al. (2001). However, the magnitude (11) should be considered only as the low limit 
for the real magnetic field of the neutron star because one need to take into account the 
contribution of the optical companion in the radiation of the accreting disk of Sco X-1 
system. Therefore the real value of the neutron star magnetic field must be increased at 
least an order and reaches the value at the level > lO*" G. 

This is the remarkable thing that three independent estimations give close values for 
the magnetic field strength in Sco X-1 system. 

It should be noted that the case of accretion disk with the chaotic (turbulent) magnetic 
fields gives the same estimation (13) as for some regular magnetic field with appropriate 
choice of the angle 9 (see formulas in Silant'ev, 2002). Besides, the model with chaotic 
magnetic fields explains naturally the " fiaring" character of circular polarization that is 
due to the fluctuations of a magnetic field. 

6. Cyg X-1: Estimation of Magnetic Field Strength in a Plasma Near the 
Black Hole 

How to extract energy from a rotating black hole is important issue in astrophysics. 
Many people have considered various alternative mechanisms for extracting energy from 
a rotating black hole. Among them the most promising one is the famous Blandford- 
Znajek mechanism (Blandford & Znajek 1977). In this mechanism a Kerr black hole is 
assumed to connect with surrounding matter with magnetic field lines. The magnetic 
field lines thread the black hole's horizon and its rotation twists the magnetic field lines 
and transports energy and angular momentum from the black holes to the accretion disk 
(Blandford 2001; Li 2002; Li & Paczynski 2000, and references therein). 

A magnetic field connecting a black hole to a disk has important effects on the balance 
and transfer of energy and angular momentum. The global structiire of black hole magne- 
tospheres involving axisymmetric magnetic field and plasma injected from an accretion 
disk has been extensively investigated for explaining varioiis observational features of 
transient X-ray binaries and active galactic nuclei (see, e.g., for the review by Beskin 



Recently Robertson & Leiter (2001) claimed an evidence for intrinsic magnetic mo- 
ments in black hole candidates. 




(15) 



1997). 
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The Electro-Magnetic Black Hole theory has been recently developed by Ruffini 
(2002). He has applied this theory to the analysis of the GRB phenomenon. He showed 
the structure of burst and afterglow of GRB can be explained within the theory based 
on the vacuum polarization process occurring in an Electro-Magnetic Black Hole. 

Unfortunately there is no till now the direct evidence of the presence magnetic field 
in near environment around a black hole. 

The marginal evidence of X-ray polarization of Cyg X-1 allows in the accretion 
disk around the compact X-ray object (black hole) (Weisskopf et al. 1977, Long et 
al. 1980). Multiple observations of Cyg X-1 were carried out by Long et al. (1980) 
with the Bragg crystal polarimeters aboard 0S08. The detection of the time-averaged 
polarization was obtained for Cyg X-1 at the level: Pi{2.6keV) = 2.4% ± 1.1% and 
Pi{5.2keV) = 5.3% ± 2.5%. If the decrease of polarization at 2.6 keV one explains as 
a result of the Faraday depolarization, the estimation of a magnetic field strength via 
Eq. (2) gives B > lO^G. It is curious that this magnetic field magnitude is quite sufficient 
to explain the famous result of measurement by Michalsky & Swedlund (1977) of variable 
optical circular polarization of Cyg X-1: Py = 5 • 10~^. 

7. Conclusion 

We estimated the magnetic field magnitudes of two X-ray binaries: Sco X-1 and Cyg 
X-1. We used for estimation two different methods. One of them is connected with cor- 
rect account of the Faraday rotation of polarization plane in the process of scattering of 
X-ray radiation (method developed by Gnedin and Silant'ev 1980,1984,1997). Another 
original method of determining the magnetic field of the objects like Sco X-1 is developed 
by Titarchuk et al. (2001). This method is based on observed quasi-periodic oscillations 
(QPO) frequencies in X-ray binaries and their interpretation as magnetoacoustic oscil- 
lations in neutron star binaries. It is remarkable that both methods give close results: 
magnetic field strength of neutron star in Sco X-1 appears at the level of Bs ~ 10^ G. 
It is wonderful that measurements of optical circular polarization of Sco X-1 give the 
same estimation of neutron star magnetic field in Sco X-1 (see Eq. (15) with taking into 
account the fact that optical radiation is generated in outer layers of the accretion disk 
comparatively far away from the neutron star itself). 

Acknowledgements. Gnedin is grateful for support from Russian Research Program: "Nonsta- 
tionary Phenomena in Astronomy" and program " Astronomy" of Russian Department of Science 
and Technology. 

References 

Beskin V. S., 1997, Phys. Uspekhi 40, 659 



10 Yu. Gnedin, N. Silant'ev & L. Titarchuk: Determination of magnetic field 



Blandford R. D., 2001, |astro-ph/0110397 

Blandford R. D., Znajek R. L., 1977, MNRAS 179, 433 

Dolginov A. Z., Gnedin Yu. N., Silant'ev N.A., 1995, In: Propagation and Polarization of Radi- 
ation in Cosmic Media, Gordon and Breach, New York 
Gnedin Yu. N., Red'kina N. P., Tarasov K. V., 1988, SvA 32, 186 
Gnedin Yu. N., Shulov O. S., 1971, Afz 7, 529 
Gnedin Yu. N., Silant'ev N. A., 1980, SvAL 6, 344 
Gnedin Yu. N., Silant'ev N. A., 1984, Ap&SS 102, 175 

Gnedin Yu. N., Silant'ev N. A., 1997, In: Basic Mechanisms of Light Polarization in Cosmic 

Media, Harwood Academic Publ., Amsterdam, 29 
Gnedin Yu. N., Silant'ev N. A., 2002, SvAL (Astron. Lett.) 28, 15 
lUing M. E., Martin P. G., 1972, ApJ 176, L113 

Kemp J. C, Wolstencroft R. D., Swedlund L. B., 1972, ApJ 173, L118 



Li L.-X., 2002, astro-ph/0202361 



Li L.-X., Paczynski B., 2000, ApJ 534, L197 
Lightman A. P., Shapiro S. L., 1975, ApJ 198, L73 

Long K. S., Chanan G. A., Ku W. H.-M., Novick R., 1979, ApJ 232, L107. 
Long K. S., Chanan G. A., Novick R., 1980, ApJ 238, 710 
Michalsky J. J., Swedlund J. B., Severny A. B., 1977, ApJ 180, L55 
Nikulin N. S., Kuvshinov V. M., Severny A. B., 1971, ApJ 170, L53 
Rees M. J., 1975, MNRAS 171, 457 
Robertson S. L., Leiter D. J., 2001, ApJ 565, 447 
Ruflini R., 2002, 



astro-ph/0202451 



Severny A. B., Kuvshinov V. M., 1975, ApJ 200, L13 
Silant'ev N. A., 2002, A&A 283, 326 

Titarchuk L. G., Bradshow C. F., Wood K. S., 2001, [a.stro-ph/0108457| 

Weisskopf M. C, Silver E. H., Kestenbaum H. L., Long K. S., Novick R., Wolff R. S., 1977, ApJ 
215, L65 



